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CORRECTION OF GEOMETRIC DISTORTIONS IN MAGNETIC RESONANCE IMAGES 

The present invention relates to methods and apparatus 
5 for the acquisition of image data and particularly, 

though not exclusively, for the acquisition of image data 
for the reconstruction of Magnetic Resonance Imaging 
(MRI) images representing an imaged object or article. 
The present invention also relates to methods and 
10 apparatus for the processing of image data so acquired 
for the purposes of image reconstruction. 

Nuclear Magnetic Resonance Imaging (MRI) enables internal 
portions of a human or animal body to be imaged non- 
15 invasively. The phenomenon of Nuclear Magnetic Resonance 
(NMR) occurs when a magnetic field B is applied to an 
object to be imaged. In the presence of the applied 
magnetic field B, the intrinsic nuclear spin of atomic 
nuclei of the material of the object are caused to 
20 precess around the direction of the magnetic field B. 



The frequency of nuclear spin precession is known as the 
Larmor frequency (o>)and is given by a> = yB where the 
constant of proportionality ( r ) is the gyromagnetic ratio 
of the nucleus and is distinctive of each atomic nucleus 
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type. Because those nuclear spins aligned with 
components of spin parallel to the direction of the 
magnetic field B reside in a lower energy state than 
those nuclei aligned with components of spin anti- 
5 parallel to B, the result is that marginally more nuclei 
will align in a direction parallel to that of B. Thus, 
the nuclei of the material of the object become partially 
spin polarised and therefore the material exhibits a net 
spin density per unit volume, p(x,y,z), pointing in the 
10 direction of the magnetic field B. 



However, the Larmor frequency (a>)of a given nucleus is 
typically shifted from the expected value o?=^ due to the 
influence of the electron cloud of the atom or molecule 
15 within which the precessing nucleus resides. The effect 
of the electron cloud is to partially shield the 
precessing nucleus from the applied magnetic field B. 
This results in a shift in the observed Larmor frequency 
(&), which shift is substantially proportional to the 
20 magnetic field B and can be used to identify the atom or 
molecule within which the nucleus resides. Consequently, 
the Larmor frequencies of nuclei within an image object 
can be used to identify and map the variations in 
molecular and chemical structures within the object, and 
25 these internal structures may then be imaged. 
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These properties are applied in NMR imaging (MRI) by 
causing the nuclei of the material of an imaged object to 
precess about an applied magnetic field B. This 
5 precession is detected, at various parts of the imaged 
object, as a varying magnetic field signal which 
oscillates at the shifted Larmor frequency (fi?)of the 
portion of the object being imaged. 

10 Additional magnetic fields known as "gradient fields" are 
then applied to the imaged object. The gradient fields 
vary linearly in the x, y and z directions of space and 
cause the precession frequency of the imaged object to 
vary from its already shifted Larmor frequency (<s>)as a 

15 function of the position of a precessing nucleus within 
the object. It can be shown that, if we denote the 
derivatives of the gradient fields in the x, y and z 
directions as G x , G y and G z respectively, then the signal 
received by the NMR device will be given by: 

20 

Signal (0 = JJJpO, y, z) • exp[z> {B + G x x + G y y + G z z)t]dxdydz = e im S{t) Eq . 1 

However, by changing variables such that: k x =yG x t\ k y =yG y t\ 
and k 3 =yG 2 t: it becomes apparent that the signal S 
25 received by the NMR device is simply given by the Fourier 
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transform of the net spin density per unit volume of the 

nuclei of the object being imaged, as follows: 

S W = Ufa** y> z> > ' GX Ptt(K x + k y y + k 2 z)]dxdydz Eq . 2 

5 

Where k = (k x9 k y9 k 2 ) is the spatial frequency vector of the 

spin density function in Fourier-Space (also known as "k- 
space") . 

10 Existing imaging techniques typically require an imaged 
object to be unvarying during or between successive image 
acquisition sequences, and this is often impossible to 
achieve in practice. 

15 This drawback is often especially problematic when 

attempting to correct image distortions arising in MRI 
images due to distortions in the magnetic field of an NMR 
device. Such images typically suffer geometric distortion 
due to several factors, such as inhomogeneities in the 

20 static field B and the susceptibility of the imaged 
object to magnetization. Such distortions arise as a 
consequence of the position of a proton in the imaged 
object being derived from its frequency of precession 
which is in turn determined according to the magnetic 

25 field in which it resides. When magnetic gradient fields 
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are applied in conjunction with the static field, a 
relationship between the magnetic field and the position 
of the precessing proton can be determined. However, if 
undesired inhomogeneities occur in the static magnetic 
field in the proximity of the precessing proton, the 
result will be a corresponding distortion or error in the 
measured position thereof, which manifests itself as a 
geometrical image distortion. 



10 As has been discussed in "A Technique for Accurate 
Magnetic Resonance Imaging in the Presence of Field 
Inhomogeneities", Hsuan Chang and J. Michael Fitzpatrick: 
IEEE Transactions on Medical Imaging, Vol. 11, No. 3, p. 
319 (September 1992), the size of the error in measured 

15 position depends upon the ratio of the magnetic field 

inhomogeneity and the magnitude of the gradient field at 
that position. 

Referring to figure 1 there is illustrated a graphical 
20 representation of a magnetic field 1 along a spatial 
direction (i.e. the x-direction) including a static 
magnetic field B component and an inhomogeneity of size 
AB. Also illustrated is a first modification 2 of this 
magnetic field arising from the addition of a gradient 
25 field G x in the x-direction, resulting in a magnetic field 
of B + AB + G x x, and a second modification 3 arising from 
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the addition of a different gradient field aG x of opposite 
sign resulting in the field B + AB + aG x x. a measured 
position x 0 derived from the magnetic field B + G x x (in 
the absence of the inhomogeneity AB) is seen to shift to 
5 position Xl by an amount: Xl - x 0 = AB/G X in the presence 
of the inhomogeneity AB, while the same position x 0 
derived from the magnetic field B + «G x x in the absence of 
the inhomogeneity AB is seen to shift to position x 2 by an 
amount: x 2 - x 0 = AB/aG x in the presence of the 

10 inhomogeneity AB. Solving these two position equations 
for the correct position x 0 is possible independently of 
the value of the inhomogeneity AB and requires only the 
two known erroneous positions (xj. and x 2 ) and the ratio 
(a) of the two gradient fields applied. This allows a 

15 correction of image distortions arising from the magnetic 
field inhomogeneity independently of information 
regarding the size or nature of that inhomogeneity. 

Present techniques implement such correction by acquiring 
20 two separate and consecutive images of an object using a 
fixed gradient field during acquisition of the first 
image, and a different fixed gradient field during 
acquisition of the second image. 



Unfortunately, NMR imaging (MRI) methods typically have 
image data acquisition procedures insufficiently rapid to 
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enable a series of images to be taken in rapid 
succession. The consequence is that movement of an 
imaged object (or changes within the object) over the 
period of time between successive image acquisition 
5 sequences of the object mean that the two resulting 
successive images are different not only because of 
magnetic field distortions but also because of changes in 
the imaged object itself. Thus, application of the above 
technique for correcting image distortions due to 

10 magnetic field inhomogeneities may fail, or be degraded, 
when applied to images acquired using existing 
acquisition methods since the two acquired images may 
well not be of the same view of the object (the subject 
may have moved or changed in between or during the two 

15 separate image acquisition intervals) . 

The present invention aims to overcome at least some of 
the deficiencies inherent in existing image acquisition 
and processing techniques so as to provide efficient 

20 image data acquisition which is less susceptible to 

movement of, or changes, within the imaged object when 
acquiring a pair of images of the object. The invention 
also aims to provide improved processing of such images 
prior to distortion correction thereof thereby to enhance 

25 the corrective effect of, or render more efficient the 

application of, the aforementioned distortion correction 
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methods and the like. it is to be noted that 
distortions in an acquired image arising from the 
magnetic field of the NMR device may be at least 
partially corrected without recourse to any pre-defined 
5 information concerning the distorting effects of that 
magnetic field according to the present invention. 

At its most general the invention proposed is to 
concurrently acquire two sets of image data during the 

10 same image data acquisition interval in which the 

magnetic field associated with acquisition of the first 
data set is caused to differ from the magnetic field 
associated with acquisition of the second data set. The 
invention proposes to repeatedly change the magnetic 

15 field in a pre-determined manner during the continuous 
acquisition of image data comprising the two sets. The 
relationship between the differing magnetic fields, most 
preferably the relationship between differing gradient 
fields, may then be employed in conjunction with the 
20 first and second data sets to correct image distortions 
arising in images constructed according to the data in 
either set alone. 

Accordingly, in a first of its aspects, the present 
25 invention may provide a method of acquiring and 

processing Magnetic Resonance Image (MRI) data from 
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Nuclear Magnetic Resonance signals generated by an object 
within a magnetic field having a predetermined spatial 
gradient, for use in reconstructing an image representing 
. said object, the method comprising the steps of: 
5 acquiring a first set of first image data items 

using a first value of said predetermined spatial 
gradient for use in constructing a first image of said 
ob j ect ; 

acquiring a second set of second image data items 
10 using a second value of said predetermined spatial 

gradient which differs from said first value thereof for 
use in constructing a second image of said object wherein 
second image data items of said second set are acquired 
before acquisition of said first set is complete; 
15 generating third image data items according to 

first image data items, second image data items and the 
ratio of said different first and second values of said 
predetermined spatial gradient. 



20 In another of its aspects of the present invention may 
omit the step of generating third image data items 
according to the first image data items, the second image 
data items and the ratio of said different values of said 
predetermined spatial gradient. The predetermined 

25 spatial gradient is preferably a gradient field G applied 
by the NMR device, and is most preferably the gradient 
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field corresponding with the scan/read direction of the 
NMR device, such as a gradient field G x in the x- 
coordinate of the NMR device along the scan direction 
thereof. 

5 

Thus, the first set of first image data items and at 
least some of (preferably all of) the second image data 
items of the second set are acquired during the same 
single image data acquisition step or time interval 
10 required for complete acquisition of the first set, 

rather than during separate or discontinuous acquisition 
intervals. By acquiring image data for two different 
gradient field values during the same image acquisition 
sequence, one obviates the need to conduct two completely 
15 separate and consecutive image acquisition 

steps/intervals , i.e. one for each complete set. As a 
result, one reduces or avoids the problems that arise 
when the imaged object moves or changes between separate 
consecutive such intervals. The first image data items 
20 and the second image data items may be stored separately 
upon acquisition thereby resulting in two distinct image 
data sets each one of which is associated with one of the 
two different values of the magnetic field gradient of 
the magnetic field of the NMR device. These distinct 
25 image data sets may be subsequently processed separately 
or together as required. 
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The second image data items of said second set acquired 
before acquisition of said first set is complete are 
preferably acquired from points in ■ Fourier-Space (k- 
5 space) which coincide with those points in Fourier-Space 
from which first image data items of said first set are 
acquired. 



Thus, prior to complete acquisition of the first set, at 
10 least some (preferably all) of the second image data 
items of the second set are acquired from the same 
Fourier-Space points as those from which at least some 
(preferably all) of the first image data items of the 
first set are acquired. Of course, if desired, some of 
15 the second image data items may also be acquired from 
regions/points of Fourier-Space from which none of the 
first image data items are acquired. 



Preferably, the method of acquiring the first set of 
20 first image data items and the second set of second image 
data items includes the steps of: 

(a) acquiring first image data items from a 
selected set of points in Fourier-Space; and 

(b) acquiring second image data items from said 
25 selected set of points in Fourier-Space; and 

(c) selecting a new set of points in Fourier-Space 
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and repeating steps (a) and (b) in respect of said new 
selected set of points until acquisition of said first 
set is complete. In this way, first and second image 
data items are alternately acquired from the same given 
5 selected set of Fourier-Space coordinates before a new 
set of points is selected and the process repeated using 
the new set until first set of image data items is fully 
acquired at least. 

10 It is to be noted that step (a) could be performed before 
step (b) , or vice versa, provided that both steps (a) and 
(b) are performed before each successive step (c) is 
performed. For example, the sequence of steps could be: 
< A ) (a) (b) (c); (a) (b) (c) ; (a) (b) (c)... etc; 

15 or it could be: 

< B > < a ) (b) (c); (b) (a) (c) ; (a) (b) (c)... etc. 

Since steps (a) and (b) require different values of the 
predetermined spatial gradient of the magnetic field, 
this means that that value must be changed when 
20 progressing from step (a) to step (b) , or vice versa. 

Sequence (A) also requires such a value change after each 
step (c) whereas sequence (B) does not. 

The selected set of points in step (a) and (c) may be 
25 distributed through Fourier-Space in any chosen manner. 
A selected set of points may form a line of points in 
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Fourier-Space (which could be part of a larger given scan 
line), or several such lines, or a curve (e.g. spiral) or 
several separate curves, or may cover one or more 
selected areas/regions of Fourier-Space each being of any 
desired shape and location. 



The selected set of points in steps (a) and (c) may each 
form a respective line of points in Fourier-Space wherein 
a set of points selected in step (c) may form a line of 
points being substantially parallel to a line of points 
formed by the selected set of points employed in 
preceding steps (a) and (b) . In this way, first and 
second image data items may be acquired alternately, line 
by line, from each one of a series of parallel lines in 
Fourier Space. The selected lines are preferably side- 
by-side and thus the image data acquisition is 
substantially a raster type scan of Fourier Space. 



Each of said parallel lines may extend in only one 
20 dimension of Fourier-space. Most preferably, this 

dimension corresponds with the read-direction or scan- 
direction of the NMR apparatus along which the device 
scans the imaged object during image data acquisition. 



The third image data items are generated independently of 
information of image distortions arising from said 
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magnetic field, such that an image of said object 
generated from said third image data items has less image 
distortion arising from said magnetic field than is 
present in an image of said object generated from either 
5 the first or the second image data items alone (i.e. the 
"first image" and "second image" respectively) . Thus, 
once the first and second image data are acquired, 
corrections of image distortions arising from the 
magnetic field may follow. The third image data items 

10 may be generated to correct both image pixel position 

errors such as apparent displacement of an image feature 
within an image frame, and pixel value (magnitude or 
intensity) errors such as apparent changes in contrast or 
the like of given image features, as between the first 

15 and second images . 

Such corrections are preferably generated after an 
initial processing of first and second image data items 
which, as discussed in more detail below, enhances the 
20 corrective effect of the third image data items. 

The third image data items preferably include pixel 
values of an image representing said object and some or 
all are generated such that the position (x 3 ) of an image 
25 pixel in an image constructed according to the third data 
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items is related to the positions (x, and x 2 ) of an image 
pixel in an image constructed according to the first and 
second data items respectively via the equation: 



oar, -x 2 

*3 ~~ 



a -I 

where a is the value of said ratio of said different 
values of said predetermined spatial gradient. 



Most preferably, the pixel intensity values ( i 3 ) of an 
image pixel at position x 3 in an image constructed 
10 according to the third data items is related to the pixel 
intensity values ( i, and i 2 \ of an image pixel in an image 
constructed according to the first and second data items 
respectively at positions x x and x 2 via the equation: 

15 Thus, the third image data items may be generated using 
only the first and second image data items and the ratio 
a , but independently of information of image distortions 
arising from said magnetic field. 



20 The value of the ratio a of said different values of said 
predetermined spatial gradient is preferably a constant 
value which may be positive (other than +1) but is 
preferably negative. The value of a is preferably -1 
(minus one) such that image acquisition requires only 
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reversal of the sign of the predetermined gradient field 
and does not require a change in its magnitude as between 
acquisition of first and second image data items. 

5 As discussed previously, pixel positions of the third 
image data items may be determined from the pixel 
positions of corresponding pixels (i.e. pixels 
corresponding to the same object feature) in the images 
generated according to he first and the second image data 

10 items. An initial processing of first and second image 
data items typically enhances the corrective effect of 
third image data items subsequently generated therefrom 
if one is able to accurately ascertain which pixels of 
the first and second images are corresponding pixels. 

15 This may be achieved by generating a "map" between 

corresponding pixels of the first image and second image 
which explicitly maps each pixel of the first image to 
its corresponding pixel in the second image. 

20 However, an initial processing step is preferably that of 
identifying and delineating the feature of interest in 
both the first and the second images.' Since such a 
feature may well be shifted as between first and second 
images, so too will the boundary/edge of that feature 

25 within the image frame of the first image relative to the 
image frame of the second image. Thus, a third boundary 
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is preferably derived according to the present invention 
in accordance with: the boundary of the feature of 
interest within the first image relative to the image 
frame thereof; and, the boundary of the feature of 
5 interest within the second image relative to the image 
frame thereof; such that the third boundary delineates 
the feature of interest when applied to the second image. 



The invention, in another aspect or in its first aspect, 
10 may include the steps of: 

generating first real-space image data items (e.g. 
pixels) from first image data items, and second real- 
space image data items (e.g. pixels) from second image 
data items; 

15 defining (using said first real-space image data 

items) a first image boundary corresponding to the 
periphery of a feature within the image frame of an image 
constructed according to said first real-space image data 
items; 

-° defining (using said second real-space image data 

items) a second image boundary corresponding to the 
periphery of said feature within the image frame of an 
image constructed according to said second real-space 
image data items; 

*5 defining a third image boundary according to said 

first boundary and said second boundary; 
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segmenting said first real-space image data items 
according to said first image boundary; and, 

segmenting said second real-space image data items 
according to said third image boundary. 

5 

In this way a feature of interest within the image frame 
of both images may be delineated or segmented from the 
rest of the image data within both images using a 
suitable image boundary determined from the first and 
10 second image data items. Since corresponding image 
features within the first and second images will 
typically shift to some degree due to the differing 
magnetic field gradients employed during image data 
acquisition, the corresponding boundaries of the 
15 feature (s) of interest within the two images will also 
shift. The third image boundary thus defines a third 
boundary which better resembles the first image boundary 
than does the second image boundary, but is shifted by a 
similar amount as the second image boundary relative to 
20 the first image boundary. Use of boundaries simplifies 
subsequent analysis of the first and second images. 

The first and the second real-space image data items are 
preferably each segmented such that: data items outside 
25 said first and third image boundary respectively are 
discarded; and, data items inside the first and third 



WO 2004/088346 PCT/GB2004/001461 

19 

image boundary respectively are retained. This removes 
noise or irrelevant data from both the first and second 
image the presence of which can otherwise seriously 
hamper subsequent image processing of the features of 
5 interest within the images. 



Preferably, the third image boundary is defined according 
to an average of the difference between the first image 
boundary and the second image boundary. This provides a 
10 smoothing of the third boundary. 

For example, the third image boundary may be defined by 
one of the first image boundary and the second image 
boundary modified according to an average of the 

15 difference between the first image boundary and the 

second image boundary. The first image boundary may be 
defined by a first image boundary vector and the second 
image boundary may be defined by a second image boundary 
vector and said difference between said first image 

20 boundary and said second image boundary may be a 

difference vector being the difference between said first 
image boundary vector and said second image boundary 
vector . 
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Each such vector being a vector preferably defined in 
image coordinates, such as pixel locations, which are 
common to both the first and the second image. 

5 Preferably, the third image boundary is defined by a 
third image boundary vector being either first image 
boundary vector or the second image boundary vector to 
which is added an averaged-dif f erence vector being an 
average of said difference vector. The value of each 
10 element of said averaged-dif f erence vector is preferably 
determined as a weighted average of the values of: 

a corresponding element of said difference vector; 

and, 

a predetermined number of elements of said 
15 difference vector which neighbour said corresponding 
element . 



For example, given a first image boundary vector f , and a 
second image boundary vector r , the difference vector is 
20 d=f-r. The averaged-dif f erence vector may be a 

weighted average wherein the j th element of the averaged- 
difference vector is a weighted average of n elements 
neighbouring the j th element such as: 
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~— s-n/2+i 



Where the weights w ± may be any suitable values including 
all w± = 1 (one), in which case the average would be a 
5 "Boxcar average" as is well known in the art. ' 

Since features of the first image will be shifted to some 
degree relative to corresponding features of the second 
image, it is typically the case that corresponding 

10 features will not be located at identical image pixel 
positions. Thus, given a third image boundary, the 
present invention may include the steps of identifying 
features of the second image within the third image 
boundary which correspond with features of the first 

15 image within the first image boundary. 

An initial processing of first and second image data 
items prior to generation of third image data items may 
preferably be performed as follows so as to map pixels of 
20 the first image to corresponding pixels of the second 
image as discussed above. The following methods of map 
generation are performed preferably only in respect of 
pixels within the aforementioned third image boundary. A 
mapping vector m is derived such that the position in 
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the first image corresponds to the pixel x 2 in the second 
image via the relation m xl :=x 2 . 



Generation of a map is performed by finding similar 
5 features between two compared parts of the first and 
second images (e.g. two compared scan profiles) . The 
present invention proposes an iterative procedure in 
which: 

Initially a mapping vector is created which maps 
10 every pixel in the first image on to the same position ii 
the second image: rn t =i; 

A measure of- similarity as between the second image 
and the mapped first image is then generated; 

Small (e.g. random) changes are made to the mapping 
15 vector (e.g. m\ =w. +1 for all i within a certain range) 
and a new measure of similarity as between the second 
image and the newly mapped first image is then generated; 

If this new measure represents an improvement over 
the previously found measure then the old mapping vector 
20 is replaced by the new mapping vector, else it is 
rejected, and the process repeats until no further 
improvements can be made. 



Accordingly, the invention in a further aspect or in its 
25 first aspect may include the steps of: 
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generating first real-space image data items (e.g. 
pixels) from first image data items, and second real- 
space image data items (e.g. pixels) from second image 
data items; 

5 (a) comparing first real-space data items with 

second real-space data items; 

(b) estimating whether or not data items so compared 
correspond with the same feature of the imaged object; 
and if not, 

10 ( c ) iteratively repeating steps (a) and (b) in 

respect of real-space data items at least one of which 
differs from any of those compared in the previous 
iteration of step (a) . 

15 This comparison may be performed collectively in respect 
of groups of first and second image pixels in which 
different groups of first and/of second image pixels are 
iteratively compared so as to identify common image 
features of varying size and/or location. The groups of 

20 the first and/or second image pixels are preferably one 
or more lines of pixels within a respective image. 
Comparison may thus be done scan-line by scan-line for 
example . 



25 The invention, in its first aspect, may comprise the 
steps of: 
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generating first real-space image data items from 
said first image data items, and second real-space image 
data items from said second image data items; 

(i) selecting a first set of data items from said 
5 first real-space image data items; 

(ii) selecting a second set of data items from said 
second real-space image data items; 

(iii) comparing data items from said first set with 
data items from said second set; 

10 ( iv ) defining a measure of similarity between data 

items so compared; 

(v) estimating according to said measure of 
similarity whether or not said real-space data items so 
compared correspond with the same feature of the imaged 

15 object; and if not, 

(vi) iteratively repeating steps (ii) to (v) in 
which at least one data item of said first set is 
compared to a data item of said second set with which it 
was not compared in the previous iteration of step (iii) . 

20 

Accordingly, each of the first and second sets of data 
items in which the pixels of the first group are mapped 
onto pixels of the second group, but the nature of the 
map is iteratively varied to find the optimal map. 
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Preferably, the compared sets of data items represent 
selected areas of a respective image containing a feature 
of interest for comparison. 

5 Accordingly, the first set and said second set of data 
items comprises image pixel values corresponding to a 
group of pixels positioned within an image constructed 
according to said first and said second real-space image 
data items respectively whereby the pixel position of 
10 each element of the group neighbours the pixel position 
of another element of the group. 

The map generation procedure may include collectively 
shifting only some pixel positions within a group of 

15 pixels upon successive iterations. This may occur when 
it appears that some of the pixels within the group (e.g. 
a sub area of pixels) have been adequately mapped, but 
others (e.g. another sub area of pixels) have not. These 
others may be separately shifted iteratively (i.e. the 

20 mapping vector changed only in respect of those pixel 

positions) in order to obtain an improved overall mapping 
vector. 



Accordingly, step (ii) of the above map generation 
25 procedure preferably includes selecting a subset of one 
or more data items from within said second set defining a 
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subgroup of pixel values whereby the pixel position of 
each element of the subgroup neighbours the pixel 
position of another element of the subgroup when the 
subgroup comprises a plurality of elements, and in which 
5 step (iii) includes comparing pixel values of the first 
set with pixel values of the second set in which the 
pixel position of at least one element of said subset is 
displaced relative to the pixel position of the same 
element during the previous iteration of step (iii) . 

10 

This may also occur simply to determine which parts of 
the mapping vector are most in need of modification. 
That is to say, by iteratively changing different parts 
of the mapping vector one can assess many different maps 
15 and determine, according to said measure of similarity 
associated with each such different mapping vector, how 
best to modify the mapping vector. 

When the estimation according to step (v) is affirmative 
20 steps (ii) to (v) and (vi) are performed in respect only 
of elements within the subset in respect of which the 
affirmative estimation according to step (v) was 
obtained. " 

25 For example, each successive subset may comprise a 

predetermined proportion of the total number of elements 
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of the previous subset (or set) from which they are 
selected. 

For example, the predetermined proportion is one half 
5 such that map generation would initially shift one half 
of the initially estimated mapping vector. The selected 
half could be shifted in one mapping direction, then in 
the opposite mapping direction, and a comparison of the 
measures of similarity obtained in respect of each such 

10 shifted mapping vector could then be made. The shifted 
mapping vector having a better measure of similarity 
could be further processed by selecting a subset of one 
half of the elements from within the previously selected 
half (i.e. one quarter selected) and repeating the 

15 shifting and similarity-comparing steps above until no 
further halving can occur. 

The measure of similarity may be a measure of correlation 
between a first set or group of pixels from the first 
20 image and a second group or set of pixels from the second 
image to which the first set is mapped via the mapping 
vector. For example, given a mapping vector in, the 
correlation may be defined according to the pixel values 
of the first set of pixels i xr and the pixel values of the 
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second set of pixels i 2 (m) chosen according to the mapping 
vector, as: 

Similarity Measure = corr(i x ,l 2 {rh)) 

5 

The larger the value of this measure, the greater the 
degree of similarity. 

Alternatively, the similarity measure may be a measure of 
10 the "Mutual Information" (MI) contained within the first 
set and the second set of data according to the mapping 
vector. This is achieved by obtaining a transformation 
m which maps the first image to the second mage such the 
information contained in both, at the same time, is 
15 maximal. The measure of similarity may be defined in 

terms of the entropy H(i,) contained in the first set (or 
subset) of pixel values, and the entropy H(i 2 (m)) contained 
in the second set (or subset) of pixel values after 
applying the mapping. Preferably, the measure of 
20 similarity is defined as: 

Similarity Measure = H(i 1 ) + H(i 2 {m))~H(i l J 2 {fn)) 

where the third term defines the joint entropy of the 
25 first and second sets. Definitions of entropy may be 
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employed such as would be readily apparent to the skilled 
person. 



Preferably, the measure of similarity is determined in 
5 respect of a given first set of pixels - such as a 

line/profile of pixels in the first image - relative to a 
plurality (e.g. three, five or more) of neighbouring 
second sets of pixels - such as a plurality of 
neighbouring lines or profiles in the second image. This 
10 helps to stabilize the map generation algorithm due to 
the averaging of noise. 



The map generation algorithm may modify the mapping 
iteratively by exhaustively searching all possible maps 

15 similar to, or in the proximity of, the given current 
map. Alternatively, the modification may be random, 
psuedo-random or directed in the direction of improving 
similarity measure. In' the latter directed variant an 
optimisation algorithm may be employed to optimise the 

20 similarity measure during said mapping modifications. 

Any suitable optimisation algorithm may be employed such 
as would be readily apparent to the skilled person e.g. 
using a Newton-Raphson optimisation algorithm or the 
like. 



25 



WO 2004/088346 PCT/GB2004/001461 

30 

The present invention may provide a method of 
constructing nuclear Magnetic Resonance images (MRI) or 
other images using the method of the invention in any 
aforementioned aspect. 

5 

The invention may provide apparatus for constructing 
nuclear Magnetic Resonance Imaging (MRI ) images according 
to the invention in its first aspect, and including none, 
some or all of the preferably features, alternative 
10 features and modified steps discussed above. 

Accordingly, in a second of its aspects, the present 
invention may provide apparatus for acguiring and 
processing Magnetic Resonance Image data as defined in 
15 Claim 24. The invention in its second aspect may be 
provided with the additional features according to any 
one or more of Claims 25 to 4 6. 



In a third of its aspects, the present invention may 
provide a computer system for use according to the method 
of the present invention in its first aspect including 
none, some or all variants/ modifications discussed 
above . 
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The present invention may also provide for the use of a 
computer system according to the method of the present 
invention in its first aspect. 



5 In a fourth of its aspects, the present invention may 

provide a program for a computer comprising computer code 
which when executed on a computer system implements a 
method according to the invention in its first aspect 
including none, some or all of the variants or 
10 modifications thereto discussed above using acquired 
image data. 



The invention may provide a computer program product 
storing a program for a computer according the fourth 
15 aspect of the invention. 



The present invention may provide an image generated 
according to, or using, and of the methods or apparatus 
of the present invention in any of its aspects. 

Non-limiting examples of the present invention follow i] 
the embodiments described below with reference to the 
accompanying drawings in which: 

Figure 1 schematically illustrates an example of 
image pixel position errors arising from distortions in 
the magnetic field of an NMR device; 
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Figure 2 schematically illustrates an imaged slice 
of an object in the coordinate system of an NMR device 
during the acquisition of image data for Magnetic 
Resonance Imaging (MRI); 
5 Figure 3 schematically illustrates alternate 

acquisition of first and second image data items from 
selected lines in Fourier-space; 

Figure 4 schematically represents a first image and 
a second image shifted relative to the first within a 
10 common image frame; 

Figures 5a and 5b schematically illustrate shifted 
boundaries of a common image feature within a first and 
second image respectively; 

Figure 6a illustrates the position of the boundaries 
15 illustrated in figures 5a and 5b as a function of pixel 
line number; 

Figure 6b illustrates the difference in boundary 
positions illustrated in figure 6a as a function of pixel 
line number; 

20 Figure 7 illustrates a map between two profiles. 

Figure 2 illustrates schematically an imaged object 
within the coordinate system of an NMR imaging device. 
The object is bathed in a magnetic field comprising a 
25 static component B aligned along the z-axis and gradient 
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fields with gradients G x , G y , and G 2 in the x, y" and z- 
directions respectively. 



The object is imaged in slices lying within the x-y plane 
5 of the coordinate system, each image slice comprising a 
raster array of neighbouring image lines ( N> prof iles") 
each formed from a line of image pixels in the x- 
direction and collectively forming a two-dimensional 
array of image pixels forming a two-dimensional image of 
10 the imaged slice. Three-dimensional images (volumes) may 
be constructed from a plurality of neighbouring 2- 
dimensional image slices. 



The imaged object may be sliced at any orientation. This 
15 is achieved by using a suitable combination of all three 
gradients Gx, Gy, and Gz. The xyz system refers to a 
reference of frame that is fixed with respect to the 
magnet/scanner of the NMR system. The user/operator of 
the system may choose to image along a plane at arbitrary 
20 rotation to the xyz planes and achieves this by using a 
linear combination of non-zero gradient fields, e.g. 
equal strength of Gx and Gy results in an imaging plane 
as drawn in Fig. 2 but with profiles running diagonally. 
In such cases a second "RPS" coordinate system is defined 
25 by: the read direction (R) ; the phase encode direction 

(P) ; and, the slice select direction (S) (2D alias multi- 



WO 2004/088346 PCT/GB2004/001461 

34 

slice acquisition) or Second' phase encode direction (3D 
acquisition) . It is to be understood that the invention 
in any of its aspects may be implemented using any 
suitable coordinate system (such as the "RPS" system) as 
5 would be readily apparent to the skilled person. 

However, in the following non-limiting embodiments for 
the purposes if illustrative simplicity, an un-rotated 
xyz coordinate system is used in which the aforesaid 
predetermined gradient field is applied along the x- 
10 dimension of the coordinate system. 

Referring to Figure 3, image data acquisition occurs by 
collecting NMR signals S(k), where: 

1 5 S(k) = jjjp(x, y, z) ■ exp[i(k x x + k y y + k 2 z)]dxdydz , 



(i.e. Eq.2 above) in Fourier-space (k-space) representing 
the spatial frequency components of the image (i.e. 

p(x,y,z)) of the imaged slice of the object. The image of 
the slice is derived from the inverse Fourier transform 
of the Fourier-space data S(k) acquired for that slice. 

The Fourier-space image data items are acquired at 
successive points within the Fourier-space over a given 
range of k x - values forming a line (scan line) in the k x - 



ine 
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dimension of Fourier-space, representing spatial 
frequency along the x-dimension in real-space, each 1 
having a different one of several fixed k y -values (e.g 
*f .*?\ *< 3) ...etc). 



For example, a scan line (e.g. Jfc« ) is selected, then: 



(a) first image data items are acquired along the 

selected scan line with the imaged object in the 
presence of a gradient field having a gradient 
value of +G X in the x-direction. Acquired first 
image data items are stored in a first memory 
store; then, once the scanning of the selected 
scan line {kf) is completed, 
(b) the gradient field is inverted to have a value of 
-G x in the x-direction (i.e. reversed relative to 
the its direction in respect of the previous 
scanning operation) . Second image data items are 
acquired by scanning the same scan line (k™) 
scanned once more with the imaged object in the 
presence of the inverted gradient field -G x . 
Acquired second image data items are stored in a 
second memory store; then, once the scanning of 
the selected scan line {kf) is completed, 
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(c) A new scan line (e.g. kf ) ) is then selected, 

having a different k value in the y-dimension of 
Fourier-Space slightly displaced from the previous 
value (k™), and steps (a) and (b) are repeated 
5 until all scan lines have been scanned. 

This process of acquiring scan lines, in which the sign 
of the gradient field is successively reversed 
alternately in respect of each of the selected scan 
10 lines, is repeated a predetermined number of times so as 
to acquire Fourier-space image data items in a raster- 
scan fashion during a single continuous image data 
acquisition time interval. Referring to Figure 3, 
"forward" scan lines (crosses "x" in Figure 3) in which 
15 the gradient field had a gradient value of +G X in the x- 
direction, and "reverse" scan lines (dots "•" in Figure 
3) in which the gradient field had a gradient value of -G x 
in the x-direction (i.e. a reversed direction), are 
acquired alternately in time such that a reverse and a 
20 forward scan is made of each selected scan line before 
the next selected scan line is scanned. Forward scan 
lines are scanned at scan periods beginning at times: ti, 
t3, t 5 , t 7 ,... etc; while reverse scan lines are scanned at 
scan times beginning at times: t 2 , t 4 , t 6 , t 8 ,... etc. Thus, 
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forward scanning and reverse scanning of scan lines is 
interleaved in time. 

The image data items S f (k) acquired in the "forward" scan 
5 lines may be labelled as such and stored separately from 
the image data items S r (k) acquired in the "reverse" scan 
lines. 

Data acquired in the "forward" scan lines ("x") alone may 
10 then be used to generate "forward" images of the object 
while data acquired in the "reverse" scan lines ("•") 
alone may be used to generate "reverse" images of the 
object . 

15 Referring to Figure 4 there is schematically illustrated 
overlapping real-space images (i.e. p(x 9 y,z)) of the imaged 
slice constructed according to the forward scan lines and 
the reverse scan lines with the image frames 6 of each 
coincident. This comparison of the forward and reverse 

20 images schematically illustrates the relative 

displacement of the position of the main feature 8 of the 
forward image of object relative to the position of the 
same feature 7 in the reverse image. 

25 The boundary 10 of the main feature 7 in the reverse 

image is defined by a reverse image boundary vector r,and 



WO 2004/088346 PCT/GB2004/001461 

38 

the boundary 11 of the main feature 8 in the forward 
image is defined by a forward image boundary vector /. 
The boundary 10 of the main feature in the reverse image 
is displaced relative to the corresponding boundary 11 in 
5 the forward image and so the reverse and forward image 
boundary vectors differ by an amount defined by the 
difference vector: d=f-r which represents the relative 
boundary displacements as between forward and reverse 
boundary vectors for the image feature of interest. It 
10 is to be noted that the aforementioned boundary 

displacement will typically be more complex than a simple 
wholesale shift of between otherwise identical forward 
and reverse boundaries, but may well include 
distortions/differences in the very shape of the boundary 
15 itself. This information will be contained within the 
difference vector d . 



Furthermore, image sub-features located within the 
boundary 11 of the forward image feature 8 will typically 

20 differ in their position - measured relative to the 
forward boundary - as compared to the position of the 
corresponding sub- features in the reverse image - 
measured relative to the reverse boundary. Thus, in 
order to correct for image distortions within the 

25 boundary of the main feature of interest, a mapping 



10 
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vector is generated as discussed below which maps each 
pixel of the feature of interest in the forward image 
onto the corresponding pixel associated with the same 
image feature in the reverse image. Once each image pixel 
5 of interest of the forward image has been paired-off with 
its corresponding equivalent pixel of interest in the 
reverse image, via the mapping vector, a third pixel 
corrected for position and pixel value can be generated 
as discussed below. 

However, initial processing of the forward and reverse 
image pixels is preferably performed to render the 
mapping and final correction procedures more efficient 
and less error prone as follows. 



15 



Initially, the boundary of the feature of interest is 
identified in both the forward and reverse images. 
Referring to Figure 4, these boundaries are 11 and 10 
respectively and are defined according to the boundary 
20 vectors /and r respectively. These boundaries may be 
identified using methods such as any readily apparent to 
the skilled person. 

Once each boundary is identified and defined, the 
25 difference boundary vector d=f-F is determined a third 
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image boundary 13 is then determined according to the 
forward and reverse boundaries, 11 and 10. The forward 
and reverse images are then segmented according to the 
first image boundary 11 and the third image boundary 13 
5 respectively. 

In this way a feature of interest (7 or 8) within the 
image frame of both images may be delineated or segmented 
from the rest of the image data within both images using 
10 the first and third image boundary. Use of these 

boundaries simplifies subsequent distortion correction of 
the images. 

The forward and reverse images are preferably segmented 
15 such that: image pixels outside said first and third 

image boundary respectively are set to a value of zero; 
and, image pixels inside the first and third image 
boundary respectively are retained. This removes noise 
or irrelevant data from both the first and second image 
20 the presence of which can otherwise seriously hamper 

subsequent image processing of the features of interest 
within the images. 

The gap 12 between a locally "inner" forward (or reverse) 
25 image boundary and the locally "outer" reverse (or 
forward) boundary, where the two boundaries do not 
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coincide, is also reduced when segmentation according to 
the third image boundary 13 is used. Furthermore, the 
third image boundary is defined according to either the 
first or the second image boundary together with an 
5 average of the difference between the first image 

boundary and the second image boundary. By using the 
average of the difference of the forward and reverse 
boundaries (11 and 10), one is able to adjust the third 
image boundary 13 to "smooth" the third image boundary as 
10 follows. 

The third image boundary is defined by a third image 

boundary vector r' being the sum (F' = F + J') of the reverse 

boundary vector f and an averaged-dif f erence vector d' 

15 which is a weighted average of the difference vector d 
such that: 

n I=0 

20 wherein the j th element d) of the averaged-dif f erence 

vector is a weighted average of the j th element, and the n 
elements neighbouring the j th element, of the difference 
vector d where the weights w A may be any suitable values 
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including all w± = 1 (one), in which case the average 
would be a "Boxcar average" as is well known in the art. 



Figures 5a and 5b separately illustrate a small region of 
5 the reverse and forward image boundaries, respectively, 
shown in the inset area 9 of Figure 4. The inset portion 
9 of each of the forward and reverse images is 
illustrated together with image pixel coordinates. As 
can be seen, due to the raster-scan nature of the image 
10 pixels in the forward and reverse images, the boundaries 
are pixelated and stepped rather than smooth. The 
forward image boundary 11 is displaced diagonally by 
three pixels relative to the reverse image boundary 10. 

15 Each Figure 6a illustrates the boundary positions of both 
the forward (+) and reverse (x) boundaries as a function 
of line number within the inset region 9. Figure 6b 
illustrates the value (circles) of the difference vector ' 
d = f-r in respect of the forward and reverse boundary 

20 vectors within the inset region 9. As can be seen, the 
value of the difference vector oscillates between the 
values of 2 and 3 pixel positions due to the stepped 
nature of the boundaries in question. The averaged- 
difference vector J' is also plotted in Figure 6b (solid 

25 horizontal line) to illustrate its smoothing effect upon 
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the difference vector d = f-r and therefore upon the 
third image boundary r' = r + d'. 

Any boundaries of relevant features of interest within 
5 any acquired image (I), whether it is within the forward 
image (I f ) or within the reverse image (I r ) , may initially 
be identified as follows: 



(1) Initially apply a Sobel filtering: Apply a Sobel 
10 kernel G** e ' separately in the x-direction (Gf Ae/ ) and the 
y-direction (Gf 6e/ )of the image I, where, 



15 



-•Sobel _ 



-l o r 

-2 0 2 



and G s , aM = 



0 0 



0 



V 



-1 -2 -1 



The final Sobel-f iltered image I Sobel being obtained from, 
I™«=I®G? M and I?*'=I®G s ; bd as a magnitude image, 



(2) then threshold the magnitude image l Sobel at a suitably 
chosen level "Th", to generate therefrom a threshold 
image I Th [i,j] such that: I Th [i,j]=l if I Sobel [i, j ] >Th; 

20 I Th [i,j]=0 otherwise; 

(3) then define "connected regions" each being the group 
of pixels in which each pixel in the group is immediately 
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next to another pixel in the group, the neighbouring 
pixel being either above, below, to the left or to the 
right. Identify the biggest "connected region" being the 
connected region having the largest number of pixels, and 
5 discard all other pixels such that I Th -* x Bigge8t ; 

(4) then "close" the image lowest by a succession of 
"dilation" and "erosion" operations: 

/aOTerf =Q B,sses ' ®K)®K where the operator 0 is the dilation 
10 operator, and the operator ® is a convolution operator. 
As will be appreciated by the person skilled in the art, 
if A®K signifies a dilation of image A with structuring 
element K, the dilation is defined as: 
C = A@K = \J(A) k 



15 where (A) k represents the translation of A by k. Thus, 
for each non-zero element k ifj of K, the image A is 
translated by i pixels in the x-direction and by j pixels 
in the y-direction, and summed onto the quantity C using 
the "OR" operator. Closing of the image may be done 

20 using a disk structuring element, such as: 



which smoothes contours, fuses narrow breaks and long 
thin gulfs, eliminates small holes, and fills gaps in 
contours; 
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(5) finally, extract the position of the boundary in line 
j by searching for the pixel position p of the first 
pixel with a value of 1 (one) in I closed [ * , j ] . Th i s i ea ds 
to the forward image boundary vector / and the reverse 
image boundary vector r e.g. such that / y = p with p being 

the position of the boundary in the line j of the forward 
image . 

Given a third image boundary, the present invention may 
include the steps of identifying within that image 
boundary corresponding/matching pixels of the first image 
and the second image. 

Figure 5 illustrates schematically a preferred method of 
matching pixel positions in the forward image and the 
reverse image. 



Initially, a first set of pixels is selected from the 
forward image, in this case a single profile *, of the 
image (i.e. one-dimensional scan line of pixels). A 
second set of pixels is selected from the reverse image, 
in this case a single profile i 2 of the image (i.e. one- 
dimensional scan line of pixels) . The pixels of the 
forward profile i, are compared with the pixels of the 
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reverse profile i 2 and a measure of similarity between the 
two profiles is made. This measure may be a correlation 
measure or a "Mutual Information" (MI) measure. Then, an 
estimation is made according to the measure of similarity 
5 as to whether or not the profiles so compared correspond 
with the same feature of the imaged object; and if not, 
the process is repeated in which at least one pixel of 
forward profile is compared to a pixel of the reverse 
profile with which it was not compared in the previous 
10 comparison. 

Accordingly, in each of the forward and reverse profiles, 
the pixels of the forward profile are repeatedly mapped 
onto pixels of the same reverse profile, but the nature 
15 of the map is iteratively varied to find the optimal map. 

Figure 7 illustrates that the intensities of the majority 
of the pixels in the forward and reverse profiles 
illustrated are identical and in the same place, but only 
20 a subgroup of pixels (representing the same image sub- 
feature) in the middle of each profile are relatively 
displaced. The correct mapping vector relating forward 
and reverse profiles is: 

m =[0, 1,2,3,4, 6,7,8,9 ,9,10, 11,121 
25 where the underlined elements are the displaced pixels of 
the central sub-feature of the profiles. 
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The map generation procedure includes collectively 
shifting only some pixel positions within the forward 
profile upon successive iterations of the above map- 
5 generation procedure. This may occur when it appears 
that some of the pixels within the profile (e.g. the 
pixels 0-4 and 9-12 of Figure 7) have been adequately 
mapped, but others (e.g. subgroup of pixels 5-8) have 
not. The unmapped subgroup of pixels is then separately 
10 shifted iteratively (i.e. the mapping vector changed only 
in respect of those pixel positions) in order to obtain 
an improved overall mapping vector. 

The mapping vector generation method according to the 
15 present example uses an optimisation procedure in which 
optimal mapping is found if the mapping vector maps the 
pixels of the forward profile into the pixels of the 
reverse profile such that the similarity between the 
mapped forward profile , as compared to the reverse 
20 profile, is maximal as follows: 

(1) find a starting mapping vector m = m Start such that 
every pixel in the forward profile i { maps onto the same 
pixel position in the reverse profile 7 2 (i.e. set m = i) ; 

(2) assess the quality of the mapping by determining 
25 a measure of similarity between the forward profile i x and 
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the reverse profile i 2 (fh) after mapping transformation 
(i.e. pixels chosen according to the mapping vector) 
using either a simple correlation such as: 
Similarity Measure = corrfi J 2 (m)) 

5 Or 

Similarity Measure = + H(i 2 (m))- H<J lt i 2 (m)) 

where the measure of similarity is defined in terms of 
the entropy Hfi) contained in the forward profile of 
10 pixel values, and the entropy H(i 2 (m)) contained in the 
reverse profile of pixel values after applying the 
mapping. The third term defines the joint entropy of the 
two profiles. The entropy of a variable X is defined as 
H (X) = ]£p(x)logp(x) summed over all values of x in the 

X 

15 vector X, and p(x) is the probability of a particular 
value x. The joint entropy of X and Y is 

H(X,Y) = J^p(x,y)logp(x t y) where p(x,y) is the probability of 
a particular value x in the same position as value y; 



20 



(3) modify the mapping vector by exhaustively 
searching all possible mapping vectors in the proximity 
of the current mapping vector (alternatively, the 
modification may be random, or directed in the direction 
of improved mapping quality, e.g. using a Newton-Raphson 
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optimisation algorithm or the like), return to step (2) 
to assess the quality of each map until no quality 
improvement can be achieved. 

5 With reference to Figure 7, the modifications of step (3) 
include the following mapping vectors for example: 
«o=[0, 1,2, 3, 4, 6, 7, 8, 9, 9, 10, 11,12, 13] 
"V s [1,2, 3, 4, 6, 7,7, 8, 9, 10, 11, 12, 13] 
«*--[-!, 0,1, 2, 3, 4, 6, 8, 9, 10, 11, 12, 13] 
10 «V = [0, 1,2, 3, 4, 6, 7, 9, 10, 11, 12, 13, 14] 
»»,._- [0,1, 2, 3, 4, 6, 7, 8, 8, 9, 10, 11, 12] 

The quality of -each is determined before modifying the 
vector to the next value and determining the quality of 
15 the modified vector. 

Initially all elements in a first half (or thereabouts) 
of an initial mapping vector are shifted to the left 
(i.e. reduced in position value) by one pixel position - 

20 the quality of that mapping vector is then determined. 

Subsequently, all elements in the first half of the given 
mapping vector are shifted to the right (i.e. increased 
in position value) by one pixel position - the quality of 
that mapping vector is then determined. The qualities of 

25 the two modified mapping vectors are compared and the one 
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with the better quality is chosen as the first new 
mapping vector. The process may then be performed using 
the first new mapping vector but in respect of the other 
half of that vector (which is identical to the same half 
5 of the initial mapping vector) so as to define a second 
new mapping vector. 

Of the aforementioned modified first half of the first 
new mapping vector, one half (thereabouts) thereof (i.e. 
10 one quarter overall) are shifted to the left (i.e. 

reduced in position value) by one pixel position - the 
quality of that mapping vector is then determined. 
Subsequently, all elements in the same one quarter of the 
given mapping vector are shifted to the right (i.e. 

15 increased in position value) by one pixel position - the 
quality of that mapping vector is then determined. The 
qualities of the two modified mapping vectors are 
compared and the one with the better quality is chosen as 
the third new mapping vector. The process may then be 

20 performed using the second new mapping vector but in 
respect of the other half of that vector (which is 
identical to the same half of the third new mapping 
vector) . 

25 The process repeats, successively halving the number of 
pixels being shifted, until no further halving is 
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possible. This represents the optimised mapping vector 
which maps the position ( x 2 ) of each pixel of the forward 
profile to the position ( x } ) of the corresponding 
equivalent pixel in the reverse profile. 

5 

Armed with this mapping one may then generate a third 
image profile from the forward and reverse profiles which 
is corrected for geometrical image distortions arising 
from inhomogeneities in the magnetic field of the NMR 
10 device as follows: 



The position (x 3 ) of an image pixel in the third profile 
is determined from the positions (*, and x 2 ) of an image 
pixel in the reverse and forward profiles respectively 
15 via the equation: 

_ ax l -x 2 

3 : 

a-1 

where a is the value of said ratio of the values of the 
gradient field employed in acquiring the respective 
forward and reverse image data (i.e. a=-l) 



The pixel intensity values (i 3 ) of each an image pixel at 
position x 3 in the third profile is related to the pixel 
intensity values (*, and i 2 ) of an image pixel in the 
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reverse and forward images respectively at positions 
and x 2 via the equation: 

i x -a-i 2 

5 The value a may be other than -1, such that not only is 
the direction of the gradient field reversed in respect 
of the forward and reverse images, but its magnitude in 
the x-direction is also changed. 

10 Modifications and variants to the above described 

embodiments, such as would be readily apparent to the 
skilled person, are contemplated within the scope of the 
present invention in any of its aspects. 



